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The method of Olver for bounding the error term in the asymptotic solutions of a second-order 
equation having an irregular singularity at infinity is extended to the general system of n first-order 
equations in the case when the eigenvalues of the lead coefficient matrix are distinct. Vector and 
norm bounds are given for the difference between an actual solution vector and a partial sum of a formal 
solution vector. Two cases are distinguished geometrically: In one it is possible to express the error 
vector by a single Volterra vector integral equation: in the other it is necessary to use a simultaneous 
pair of Volterra vector integral equations. Some new inequalities for integral equations are given in 
an appendix. 



1 . Introduction 

Let us first establish our notation. To refer to an m X n matrix A we shall use A = [a/j] (i= 1, 
2, . . . , ra7=l,2, . . . , n), while to refer to the (7, 7)th element of A we shall use {A}/j = a,j. 
We shall write |A| for the matrix of absolute values of the elements of A, and for A and B real and 
of the same dimensions, A ^ B will imply that {A},j ^ {B}/j for all (7,7"). All the integration paths 

dt{s) 
ds 

is continuous and non-vanishing. The notation _ | \{t) | where \(t) is an m X n matrix will be 

used to denote the mX n matrix of non-negative numbers whose elements are the least upper 
bounds on 3P of the corresponding elements of |V(l)|. 

The system of differential equations we consider takes the form 



we shall use consist of a finite number of Jordan arcs t — t(s) {a ^ s ^ b) on each of which 



dW = ,. 
dz 



S ^A k z- k + %(z)z-°\w (1.1) 



where r and s are non-negative integers. Each A* is a constant nXn matrix. The n 2 elements 
of ?ls(z) are holomorphic in some domain 2 (which may be part of a Riemann surface) extending 
to infinity. For any path 9* extending to infinity in 2), we suppose that ^ | {2i s (z)}ij | is finite for 

7,7=1, 2, . . ., n. For notational convenience we shall write %o(z) = ?I(z), A_j=0. The integer 
r denotes the rank of the system (1.1). 

It is known (|2) 3 ) that the system (1.1) has a linearly independent set of n formal vector solu- 
tions which are asymptotic expansions of actual solutions of (1.1) as z—> °° in 3. The purpose of 
the present paper is to establish bounds for the difference between the actual solutions and the 
partial sums of formal solutions. This is achieved by an extension of the error analysis developed 
by Olver [1] for second-order differential equations. 



'An invited paper based on the author's Ph.D. thesis, University of Alberta (1965). 

2l Present address: University of Michigan, Ann Arbor, Mich. 

3 Figures in brackets indicate the literature references at the end of this paper. 
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2. Transformations on the Original System and Formal Solutions 



For purposes of obtaining error bounds we require explicit rules for constructing formal solu- 
tions. In this section we therefore construct transformations which convert (1.1) into canonical 
form, and then we construct formal solutions for the transformed system. 

THEOREM 2.1: There exist constant nXn matrices To, T\, . . . , T r +i such that the transformation 



X = £(z)W = 2 T k z- k W 



(2.1) 



reduces (1.1) to the canonical form 



dW 

^Wg(z)w, 

dz 



(2.2) 



where S(z) can be expanded in the form 

£(z) = z r+1 ^ C k z- k + z- 2 2 B k z- k + 93 s (z)z- s . 



k=o 



s-l 
k=0 



(2.3) 



Here the C k and B k are constant n X n matrices, C k being diagonal and 33 s ( z ) is an n X n matrix 
whose elements are holomorphic for all sufficiently large z in 3, and uniformly bounded on a path 
3P extending to infinity in 3). 

We prove this important theorem here in order to furnish an explicit algorithm for computing 
the coefficient matrices Ta : (A = 0, 1, . . ., r+1). For an earlier proof see Birkoff [3]. 

Since A has distinct diagonal elements there exists a nonsingular matrix T such that To ^oTo 
= Co, where Co is a diagonal matrix with distinct diagonal elements. Accordingly we assume that 
the transformation X = T W has been made and that A is already in diagonal form. Then we 
may take T = I in (2.1). On making the transformation (2.1) in (2.2), we obtain 



dz 



=z r Z-\z) 



%z)^Xz)- 



d%(z) 



dz 



]w. 



We choose the matrices Ti, T 2 , 



., T r +i so that the matrix 



G(z) = £-»(*) \%{z) m -z-' ^] 



has the form (2.3). For this purpose we need only consider the coefficients of z _1 , z~ 2 , 
z" r_1 in the equation 



5t(z)©(z) = Sl(z)S(z)- 



dz 



(2.4) 



(2.5) 



(2.6) 



Expanding (2.6) and equating equal powers of z, we obtain 

2 CT,_ fc C ft - A,T,_ A .) = O (/ = 0, 1 , . . . , r + 1). 



With / = 0, (2.7) gives C = Ac With Z = 1 we have 

T,C - A T, + Ci - A, = O. 
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(2.7) 



(2.8) 



Now for arbitrary Ti, the diagonal part of T1C0 — AoTi = TiAo— A0T1 is zero. Thus Ci = diag 
(Ai). The off-diagonal elements of Ti can be determined by the solution of 

{T 1 } <{Aob-{Ao}u) = {A 1 } u 

under the assumption that the diagonal elements of A () are distinct. We set the undeterminable 
diagonal elements of Ti equal to zero. 

Assume that for some positive integer / ^ r we have solved for Ci , . . . , C/ and the off-diagonal 
elements of Ti, . . ., TV and that we have set the diagonal elements of Ti, . . ., T/ equal to zero. 
Then from (2.7) we have 

TViAo-AoT z +i+ J (Ti+i-*C*-A*Ti + i-*) + Ci + i-Ai+i"=0. (2.9) 

A = l 

Again, the diagonal part of T/ + iA — A T/ + i is zero. Thus the diagonal part of 

2 (T, +1 - fc C fc -A fc T /+1 _A:) + C/ + i-Az +1 (2.10) 

k = \ 

is also zero. This enables us to determine Cj+i. The off-diagonal elements of T/+i are determin- 
able from the scalar equation 



{Ti+i}«({Aob-{Ao}«)+ X (T/ + ,-a-C a --A,T / + 1 _,)L = (2.11) 

since the diagonal elements of A () are distinct. The diagonal elements of T*+i are again not deter- 
minable and we set these equal to zero. This proves by induction that each T k - and Ca (k = 0, 1, 
. . . , r+ 1) can be determined as stated in the theorem. 

To determine B^ and 93*(z),we again expand (2.6) and equate coefficients of z~ r ~ l ~ k , k^l. 
We obtain 

min(A:, r+1) 
BA- = AA- +r + 2 — 2 (TfBfc-j — Afc + r + 2-zTj) 

/=1 

- 2 (TlC* + r+f-I-Ajb + r+l-lT|) + *Tft (2.12) 

l = k+\ 

for all k i^ 0, where we define Tjt — for fc>r+l. Since the elements of 91 (z), S£(z), -t _1 (z) are 
holomorphic and bounded for all sufficiently large z in 2) , the elements of 93*(z) are similarly holo- 
morphic and bounded. 

This completes the proof of Theorem 2.1. 

Theorem 2.2: The system (2.2), in which ©(z) has the canonical form (2.3), possesses a formal 
independent series solution matrix of the form 4 



W(z) = U(z)expG(z) (2.13) 



4 Here and henceforth exp B = V -jt for any square matrix B. 
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where O(z) is diagonal and has the form 

£(*)= t ^±Z z r+ '- k + Qr + . lnz= 2 ^^+ Cr+1 Ir 
am/ 5 



(2.14) 



U(z) = I+2U k z- k . (2.15) 

k = l 

PROOF: Substituting (2.13) into (2.2), canceling exp <Q(z) and equating coefficients of equal powers 
of z" 1 , we find that 

min(/t,r+l) k-r-2 

X (U*-jQ^-QU*-i)-(*-r-l)Ujk.r.i - Y BjUft-r-i-j-CKife-O, 1, 2, . . .) (2.16) 

j=0 j = 

provided that we define \Jj = for ;< 0. The proof of Theorem 2.2 will be complete after we 
prove the following lemma. 

LEMMA 2.1: With U = I and G(z) diagonal, equation (2.16) uniquely determines Uj and Qj, 
7 = 0,1,2, .... 

PROOF: With /c = we have 

Qo-Co. (2.17) 

With k=\, we have 

UiCo-C Ui + Qi-Ci = 0. (2.18) 



Now if Ui is an arbitrary nX.n matrix, the diagonal part of UiC — C Ui is zero. By assumption, 
Qi and Ci are diagonal; hence Qi = Ci. Since the diagonal elements of Co are distinct, the off- 
diagonal elements of Uj given by the solution of {Ui}/j({C }jj— {C }n) = are zero. The diagonal 
elements of Ui are not determinable at this point and all we can say here is that Ui = Di, where 
Di is an arbitrary diagonal matrix. 

Assume that for some integer k in the range 1 ^ k ^ r we have found Qj=Cj and Uj=Df 
where D, is an arbitrary diagonal matrix and 1 t^j^k k. Then solving for U*+i and Qk+i we have 
from (2.16) 

Ua- + iCo-C U^i + 2(D a . +1 _ j C j -C j D a . +1 _ j ) + Q a - +1 -C. +1 =0 (2.19) 



or, since diagonal matrices commute, 

U A -fiCo-CoUA- + i + OA-fi-A A - + i-0. (2.20) 

On comparing (2.20) with (2.18) it follows by exactly similar arguments to those used for (2.18) 
that Qa + i^Ca + i and that U A +i = D A +i, where D A + i is an arbitrary diagonal matrix. We have 
thus proved by induction that for every integer k in the range ^ k ^ r-\- 1, Q A = C A and U A — D A -, 
where the D A are undetermined diagonal matrices. 



5 The symbol "=" here and henceforth denotes a formal equality. 
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With k = r + 2, we have, from (2.16) and the above results, 

U^Co-CoU^+XlDrri-jC—QD.^-^ + BoD.-O. (2.21) 

or, since the sum on the left is a sum of diagonal matrices, 

U r+2 C - CoU,. +2 - Bo - D, - O (2.22) 



The diagonal part of Ur+sCo "" C©U r + J is zero; hence Di = Ui= — diag B . The off-diagonal 
elements of U r +2 are given by the solution of {U r +2}/j({Co}jj — {Co}/») = {B }/j. The diagonal part 
D r +2 of U r +2 cannot be determined at this point. 

Assume that for some positive integer /, the matrices Ui, Us, . . . , U/ as well as the off-diagonal 
elements of the matrices U/+i, Ut+s, . . . , Ur+M-i have been determined while the diagonal mat- 
rices consisting of the diagonals of U*(/+ l^Ar^r-hZ-hl) are still undetermined. Then, from 
(2.16) we have 

Ur+H-lCo-CoUr+l+l+2 (Ur+I+I-iQ-C|Ur + l+I-|)+y B/-iU| + i-j-(/ + l)U|+i»0. (2.23) 



On writing Uf — Uj + Dj, j > /, where the diagonal elements of Uj are zero, and using the fact 
that diagonal matrices commute, we find that 



Ur+l+lCo-CoUr+l+1 + 2 (Ur+l+I-jQ-QUr+l+l-iJ-J B j -,U/ + ,-~(/+ 1)(U/ + , + D, + ,) = 0. (2.24) 

j=l j=l 



Thus each Uj, /SI, is uniquely determined. This completes the proof of Lemma 2.1 and of 
Theorem 2.2. 

Clearly, the above procedure can be modified to compute a particular formal series solution 
vector independently of the others; the y'th series solution vector is 



Wj(z) = JJ j(z)e«j (z) (2.25) 

where U/(z) is the 7th vector in U(z) and qj(z) = {d(z)}jj. For later convenience we set 

«»- t Jt k _L ' r+1 - k + <U.r+i log z. (2.26) 



The representation (2.25) illustrates that we have a one-to-one correspondence between the 
eigenvalues of A and the formal vector solutions. 
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3. Error Bounds for the Formal Partial Sum Approximation 
3.1. The Differential Equation for an Approximation 

We denote the yth column of Ujt (A; = 0,1, . . .) as defined by Theorem 2.2 by Uf*. Starting 
with the partial sum 



^WH^V-'e'^ (3.1) 



we define a vector Rj m (z) by the differential equation 

d$> jm (z) 



z r &(z)% m (z) = RjM) e *f* (3.2) 



where (£(z) is defined by (2.3). 

Expanding (3.2) and using equations (2.26) and (2.16), we obtain 



r+m+l rmin ( M , r+1) 

R Jm (z)=- £ £ (Qk-lq^UU + in-r-DUtv-r-i 

fx=o L ^=0 

HA— r— 2 m-l 

+ ^ B*UU- r _ 2 _ fc z'-"- ^(33 m _ fc ( z )UA)2- m - 2 (3.3) 

where Ufjfc = Ujjfe if ^ k^ m — 1; O otherwise. 

Consider the ith element R ? j m (z) of Rj m (z). If i=j then using (2.16) again we obtain 

m— 1 

R /JW (z) = niU jjw z-»- J - 2 {(Sm-^Uj,)}^— 2 (3.4) 

A-=0 

where {\j}i denotes the ith element of the vector Vj. 
We also define 

qij{z) = qi{z)-q j {z). (3.5) 

For i ¥"j we again make use of (2.16) in (3.3) to obtain 

Ri jm{z) = q [j{z)Vi jm z- m - z[q [j{z)- (q i0 - qj o)z r ] Ul jruZ'^ 1 

"min (^., r+1) 

2 

k=/ji+l-m 

tx-r-2 

+ I 

k= 

where Uijk is the ith element of the vector Ujk- 

3.2. Two Possible Cases 



fx=m-hl 



m+r+i min ju, r+ij 
— 2 2 (^"~^*)U|j, M -fc + (/Lt— r— l)Utf fA -r-i 

a— r— 2 m-l 

fe=Q J A:=0 



In this section we shall set up integral equations for the error vector, that is. the difference 
between an actual solution vector of (2.2) and a partial sum of the formal solution vector (2.13). 
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Our aim is to construct actual solution vectors for which the formal solution vectors obtained in 
section 2 are asymptotic expansions, as |z|— » oo. For this reason, and for purposes of obtaining 
good error bounds we assume 3f is such that we can choose our fixed end-points of integration, 
which we denote by £*, at oo. We shall show that under these conditions it is possible to express 
the error vector by at most a simultaneous pair of Volterra vector integral equations. 

The error bounds we shall obtain are considerably sharper in the case when we can express 
the error by a single Volterra vector integral equation than in the case when we require a simul- 
taneous pair of Volterra vector integral equations. It is thus desirable to use a single equation 
whenever possible. 

We have already noted at the end of section 2 that with our construction of formal solutions 
there is a one-to-one correspondence between the eigenvalues of Co and the vector solutions. 
Let \j be an eigenvalue of Co. 

DEFINITION 3.1: The eigenvalue Aj will be called an extreme eigenvalue if there exists a path 8P 
lying in Q) , joining some point z(|z|< oo) in 2 with £k(|£kl = °°) suca that Re qij(t) (see equation (3.1)) 
increases monotonically 6 as t traverses 2? from z to £ k , for i=l, 2, . . . , n. Otherwise Xj will be 
called an interior eigenvalue. 

Consider the eigenvalues of Co as points in the complex plane enclosed by the smallest pos- 
sible closed convex polygon 11. Except for rotation, the polygon 11 and the differences between 
the various eigenvalues of A are left invariant under the transformations W = \e Xz , W= 3^(z)X 
(5£(z) as in Theorem 2.1), z = co£, co a constant such that |co|=l. Thus we may assume, without 
loss of generality, that the vector solution of interest corresponds to the eigenvalue Aj = 0. More- 
over, in the case when zero is an extreme point of the polygon II, we may assume that all other 

7T 

eigenvalues of A are within the closed sector |arg A|^= a ^ — and that there is at least one eigen- 
value on each of the rays arg k = a and arg \ = — a. In the case when zero is an interior point of 
II, we assume that no eigenvalue (other than the zero eigenvalue) is located on the imaginary axis, 
that there are some eigenvalues of A in each of the sectors |arg X|^ a and \tt — arg X|= a, = a 

77 

< — , and that there is at least one eigenvalue on each of the lines arg \ = a and arg X = — a. 

Suppose for example that we have the case indicated in figure 1, that (^(z) is regular for all 
|z|>p, that r = and that q,, , + i = qn=0 (i = l, 2, 3, J). It is then clear that |eV| (i=l, 2, 3,;) 



increases monotonically along any path for which 



arg 



dz 

dx 



— — a. Thus assuming the points 



A; (i=l, 2, 3,7*) to be in the z-plane as indicated in figure 1, we can, for example, connect any point 



8 For any two points z,, z 2 in the order z, z,, z 2 , £ fc on <?> we have Re gyfo) % Re gy(z 2 ). 



Figure 1. 
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3 77 

in the interior of the overlapping sector |arg z\^~ tt — a — 8, ^ § < — of figure 1 with +<*> by a 

4 4 



path 8P consisting of at most three straight lines along each of which 



arg 



dz 

dx 



— — a and none of 



which passes through the origin. 

For arbitrary rank r+1 it follows that, since |arg qij{t) — arg [(#/o — qjo)t r+1 /(r-\- 1)] | can be 
made arbitrarily small by taking \t\ sufficiently large, we can connect any z (\z\ sufficiently large) in 



argz- 



2irk 



r+1 



^(t-"- 8 )! *^!)'*-^ 1 ---^ 



to 



exp 



r+1 



by a path & along which Re q$(t) increases monotonically for i = 1, 2, . . . , n, provided that 



a< 



IT 



In the case when a = — there may still be a region between two parallel straight lines (no 

longer a sector of positive angle) in which there is at least one path & joining z with £a- along which 
Re qtj(t) increases monotonically for all i. 

Let us now consider a case in which kj is an interior point of the polygon. If the eigenvalues 
of Co are as indicated in figure 2 then, given any number 8in0^i8<7r — a there is a path joining 

77 

arg z — — ^ 77 — a — 8 and + °°, consisting of at most four straight lines 



°°, any point z in the sector 



(none of which passes through the origin) such that Re (Xi£), Re (k 2 t) increase monotonically while 
Re (\zt) decreases monotonically as t traverses S* from — °° to °°. If we require that a neighborhood 
of the origin is to be avoided by ^, \z\ may need to be taken large when 8 is small. Similarly, it is 
easy to see that there is such a path joining — °°, any point z in the sector |arg z -\- 2 tt\^ rr — a — 8 
and -h 00 . Thus these two families of paths, one of which passes around the origin in a positive 
sense and one in the negative sense as t traverses 8P from — °° to + °°, cover the complete neighbor- 
hood of infinity. 

In the general case when Xj is an interior eigenvalue of II, we divide the integers N= 1, 2, . . . , 
n into two disjoint classes Ni and N 2 such that if ki is in |arg k\ ^a then 1GN1 while if A.* is in 
\tt — arg k\ ^a then i E N 2 . Note that j may be taken either in Ni or in N 2 . Clearly there are 
two families of paths connecting ^ 1) = exp [(2 A: ± 1)tt V— l/(r + 1)], and ^ 2) =exp [2A:ttV— l/(r+l)] 




Figure 2. 
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(k — 0, 1, . . . , n) such that as £ traverses ^from £ ( A 1} to £ ( A 2) Re qij(t) decreases or increases mono- 
tonically according as i E Ni or i E N 2 . These two families of paths cover the sectors | arg z — (2k 



±1) 



r+1 



r+1 



THEOREM 3.1: Let k } be an eigenvalue ofC {) . Ifk^ is a vertex of II, then it is an extreme eigen- 
value; if it is an interior point ofU then it is an interior eigenvalue. 

Let S(z) be holomorphic in a domain 2 extending to infinity in the direction arg z = 27rA;/(r+l), 
for some integer k in ^ k =i r. Let <? be the smallest closed disc with center at the origin con- 
taining the zeros of q\j{t) (i= 1, 2, . . ., n; i #y'4, and let 2* = 2 — %— {°°}. 7 We define a region 
3}(z, £/<•) to be the union of all points z such that there is a path 3P connecting z and £&= exp [27r/c V— 1/ 
(r+1)] satisfying the following extreme eigenvalue conditions:* 

(1) Except for £*, ^ lies entirely in 2)*\ 

(2) For any two points t\ and £2 in the order £*, £1, £ 2 , 2 on ^,we have 



|exp[^ (^) — qij(ti)]\ ^ 1 (1=1,2, . . ., n); 
% (t~ 1 ) = \ \t~ 2 dt\ is bounded. 

The variation symbol introduced in condition (3) above is more generally defined as follows: 
If F = F(z) is a vector of holomorphic functions of z, we define 



nm 



)r ' )oj>\ 



dF 
dz 



dz 



Let us set 

Dj(z) = €l(z)-qj(z)l. 
If a vector Wj(z) satisfies (2.2), then by (3.2) the error vector 

e im (z)={W j (z)-^j m (z)}e-"J {z) 
satisfies 

— € Jm (z) - Dj(z)€ jm (z) = -<8(z)€ jm (z) - Rjm(z). 

where 23(z) = SBofz). Now if €j m , k(z) is a solution of 

€jm,*(2)=P e D ^^[r 2 m)€jm, k(t)-Rjm, k (t)]dt 
Jik 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



where the path of integration is chosen as described above, then €j m ,k(z) also satisfies (3.10). 

If kj is not an extreme eigenvalue, let 2l(z) be holomorphic in a domain 3 which contains, or 



extends to infinity in, one of the sectors 



arg z — (2k± 2) 



r+1 



— lor some integer k in ^ k ^ r 



(r+1) 



7 {°°} = {e w |-°o=£0<°o}. 
8 Compare [1], sec. 5. 
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and all z sufficiently large. We reorder, and if necessary relabel, the elements of the equation 
(3.10) so that we may write them in the form 



A 
dz 

A 
dz 



j z €$(z) - D^)€$(z) = z- 2 [a, (z)e$(z) + a 2 {z)e%z)'\ - R$(z) 



efAz) ~ DW(z)e$(z) = z~ 2 i[a 3 (z)e$(z) + Cu(z)e$(z)] - R$(z) 



(3.12) 



where the top line contains all rows of (3.10) such that iGNi and the bottom line contains all 
rows of (3.10) such that iEN 2 (see page 174 ). The diagonal matrix D/z) in (3.10) is accordingly 
subdivided into two diagonal matrices Dj x) (z) and Dj 2) (z); the matrix 93(z) is accordingly partitioned 
into four blocks ct/c(z), (A = l, 2, 3, 4) and Rj m (z) is split into two vectors in the same manner as 
£jm(z) was split. 

Let & be the smallest closed disk with center at the origin containing the zeros of \%v) -(i= 1, 
2, . . . , n\ i 9^j) and let 3)* =2 — % — {°o}. We next define a region @((ft\ & 2) ) to be the union 
of all points z {z^ t$\ z ^ £i 2) ) such that there is a path (P connecting Qp 9 z and £l 2) (in that order) 
where ffl is one of the points °° exp [(2k ± l>7r V—l/(r + 1)] , and £k 2)==0 ° exp [27rA: V— l/(r+l)], 
and satisfying the following interior eigenvalue conditions: 

(1) Except for #J> and £g>, ^ lies entirely in ^* ; 

(2) For any two distinct points t\ and t 2 in the order Q}\ t\, t 2 , £$\ on (P we have |exp [qij{t 2 ) 
-qij{ti)]\ ^ 1 if i E Ni, |exp [«#i)-«j(fe)]| ^ 1 if i G N 2 ; 

(3) With ai, a 2 , a3, and cm defined by equation (3.32) below, 2 {a\ + a 4 + V (ai — a,4) 2 + ^0,20,3} 

T^ (r 1 ) < 1. 

It follows that the vector of analytic functions which satisfies 



€$, k (z) 



(I 







,DS'»fe)-D5'i(«) 



I 



,(2) 




ai(0 a 2 (*) 
a 3 (*) a 4 (t)_ 



_e$,k(t)_ 






<fe, 



(3.13) 



where the integrals are taken along a path 8P as described above, is a solution of (3.12). 

3.3. Treatment of Rij m (z) 

\ 
In this section we examine the integral 

Rr jm (z)=t Z e^ z) -^R ijm (t)dt 



(3.14) 



where the path 3P of integration satisfies either of the sets of conditions described in the previous 
section and £ is a point appropriately chosen at infinity such that as t traverses 3P from £ to z, 
Re qij(t) decreases monotonically. The equations (3.4) and (3.6) show that when i=j\ R ijm (t) 
= 0(\t\- m ~ 1 ), |*|-*oo, but if iy*j, then R ijm (t) = 0(\t\ r - m - 1 ), |f|->oo. For this reason we cannot 
proceed directly as in [1] to bound the error. 

When i #7 we first integrate (3.14) by parts to obtain a good bound for Rf jm (z). The equation 
(3.6) may be written 
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R ijni (z) = ~ e <hM) j- z (e quiz) U ijm z->») - z 



q'jj(z)—(qu) — qjo)z r + 



m 



IT.. 7 -m-\ 



TO+r+1 Tmin (/x, r+l) 
^t=m+l L A'=/x+i-m 

lx-r-2 m-\ 

+ J) {B,U,, M ,.,_,,} | . *-»- ^ {(*m-fc(z)Uj fc )} f 2-" 1 " 2 . (3.15) 

The first term on the right of (3.15) is already in a form suitable for integration by parts. Consider 
now the polynomial P,>,U) consisting of all of the right-hand side, excluding the first and last terms. 
It is clear that the dominant term of P,>,(z) is of order z r ~ m ~ l as |z| — ><». On substituting Pij m {t) 
for Kijm(t) in (3.14), integrating by parts and denoting the resulting integral by P*„„ we obtain 



p* ( r \= jiggfe) i 

•>'" {Z) qlM 



J i idt\ q[j(t) ) 



(3.16) 



Collecting terms, we obtain 



Riw»)=r 



and since 



mU,„,r'«- + ( e "" , ^" u, ' , -l)|(^) -e^-^f {%,,-MVj,}^"-^ dt, 



(3.17) 



, iiijUi-iul'l \ < - , CL (Pijm(t)\ ... i 



*l-'"- 2 ), 



|*|— > oo, the right of (3.17) is readily bounded. In the next section, we use this bound in obtain- 
ing vector and norm bounds for €j,„{z). 



On the other hand for i=j\ we see from (3.4) that 

fz V TO - 1 

Rfj,„(z)= mUjjmt-™- 1 - X {%m-k(z)Uj k }jt- m - 2 
and obtaining a bound for this integrand is even easier than for (3.17). 



dt 



(3.18) 



3.4. Error Bounds for the Extreme Eigenvalue Case 

We shall first obtain a vector bound. With t\ and t 2 in the order £*, t u * 2 , z on ^, we define 
sup 



(/n+l)yym = 



ri, ^ E ^ 



^V'W 



t w+2 



rf / Puro(£) \ 



^ (Y ym(t)\ t , n + 



(m+l)y jjm - t G ^ 



m-1 

2 {»— mou^ 



(3.19) 



(»*/) 
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We also denote by yj m the vector with ith element y«»; i.e. {yjm}i = yijm, and define 



B = 



sup 



e^-^i>{» (*!)}«. 



, i, 5 = 1, 2, . . . , n. 



(3.20) 



In the part II of the present paper we shall show that there exists a unique vector €j m ,k(z) of 
holomorphic functions which tends to zero as z—*£k and satisfies (3.11). The elements of €j m ,k(z) 
are holomorphic in a domain (actually a Riemann surface) which includes the region £&{z, £a) 
defined by the extreme eigenvalue conditions in Section 3.2. 

On combining equations (3.19) and (3.20) with (3.11) we obtain 



Icfejkfr)! Si [U- s Bcjm,, jk(t)| + wt|Uj»i-" , - 1 | + (m + l>»i«t- , "-«l|€fc| . (3.21) 

Let a be the value of s for which £(s) = z t(s) as defined in the Introduction. We apply Lemmas 
1 and 2 of the appendix with E(cr)=l, ¥(s)ds=B\t- 2 dt\ 9 <p(s)= |cjto,*(£)|, G(5)^5=(^|Uj m ^ m - 1 | 
+ (m + l)y jm |j- m - 2 |)d£ to obtain 

THEOREM 3.2: Corresponding to an extreme eigenvalue of Xj o/Co, tAe equation (2.2) possesses 
an actual solution vector W jmjk (z) depending on £ k and an arbitrary positive integer m 5ac/i f/ia£ 



W Jm . k (z) = 



m — 1 

^ U js Z" s + € jm , k (z) 
L S = 



e <7j(z) 



(3.22) 



where 



|€ Jm , k (z)| ^ exp{B^ (t-^ilUjWp (t- m ) + y jm ^ (t--" 1 )}^ z E 0(z, £k) 



(3.23) 



and «^(z, £ k ) t5 defined in section 3.2. In (3.22) eac/i vector U jk (k = 0, 1, 2, . . .) is *Ae ]th column 
vector of the matrix U k defined by Theorem 2.2; q^(z) is the ]th diagonal element of the diagonal 
matrix O(z) (equation (2.13)). In (3.23) B is an n X n matrix of non-negative elements defined by 
equations (2.3) and (3.20), while the \th element (i=l, 2, . . . , n) of the vector y m is defined by 
equation (3.19). 

Let us now obtain a norm bound for €j m ,k(z) in equation (3.22). For this purpose we define a 
vector Vj(^2, ti) as follows: Let the \th element (m + l)V/j(^2, t\) of (m+ l)V/*2, t\) be the number 
inside the absolute value signs on the right o/(3.19). With t\ and t 2 defined as for equation (3.19) 
we define 9 



B- 



sup 



$ 



II830II 



^"* lf "e^ |Vi(fe,l) 



Theorem 3.3: A norm bound for the vector e jm , k (z) in Theorem 3.2 is given by 
||€ jm , k (z)|| S exp{B^ (t-»)}{||U to ||^. {t- m ) + y- m n (t— •)}. 



(3.24) 



(3.25) 



9 Here and henceforth we assume that compatible matrix and vector norms are taken, i.e., ||BV||^ ||B|| ||V||., 
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3.5. Error Bounds for the Interior Eigenvalue Case 



Put t 2 = z, ti = t, and let \j(z, t) be denned as for (3.24). We split the vector V/(z, t) into two 
vectors \ { j l \z, t) and \f\z, t) in a manner corresponding to the way in which Rj W! (0 was split in 
(3.13). We analogously split the vector U,,„ into U$ and U^, such that the equation (3.13) may be 
written 



e jm,fA z ) 




h> ] 





i 


r D/»(2)-I>/'(/) ft 

e 1 


1 


;~ 2 


.«&*(*)- 




t ° /iJ 


i 


e j 




"ai(0 


a#)~ 


r«KU*>i 






X 








4- 










_n :{ (/) 


a 4 (0_ 


U 2) j(*)J 

L jm,lc '- 1 




- mU$t 



' m - l + (m+l)\ljKz 9 t)t- m ~ 2 



BllJ^-' R - 1 + (/Jl+ 1)V ( ; 2) (2. /)/ 



A. (3.26) 



Again we shall first obtain a vector bound. To achieve this, we define 



a*= f s e u yiw<)} is i] 



(4-1,2,3,4) 



CtfW^mlUjftl^-H + 'i.fce^m'+^VjJil^fOlk 



(3.27) 



where in the last of (3.27) ti and t 2 are on ^ in the order £i, *i, £ 2 , $2 for s= 1, and {1, f 2 , £1, £2 for 
5 = 2 respectively. 

DW*)~*DWf) 

In addition, each element of the diagonal matrices e ' ' (5= 1, 2) on the right of (3.23) has 
an upper bound of 1. 

In part II of the present paper we shall show that there exists a unique nXl vector of functions 
satisfying (3.26) such that c^ A .(z)— >0 as z — > $P, e^j A .(z)— >0 as z^> ^ 2) . Furthermore, each ele- 
ment of this vector is holomorphic on a Riemann surface which includes the region &{Ql\ Qr } ) 
defined by the interior eigenvalue conditions of section 3.2. 

On substituting the above results in (3.26), we obtain 



Ll«$»l. 



/" 



ill il J 



03 »4_ 



«$>)l" 


|r»| + 


"C ( „V(0 


e$»IJ 




Lo|)(t). 



Iflfel 



(3.28) 



On transforming to real variables as for (3.21), using Lemma 5 of the appendix and transform- 
ing back to complex variable notation, we obtain 

THEOREM 3.3: If, for an interior eigenvalue Xj of C () , the interior eigenvalue conditions of 
section 3.2 are satisfied, then the equation (2.2) possesses an actual solution vector 



rnut) 

lw^, k (z)j 



m-7 


Wl 




r<»] 


2 




z" s + 




s = 


W 




Us»J 



#> 



(3.29) 



|< k (z)|g exp {a 1 ^», t (t- 1 )}«^ ) (4 1) . i) + C<(fl» z)[exp {a^ (t"*)}^^, ffO + B], (z6^»,& 2) )) 
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and 

#2(ii, v)= P Cg(t)|dt| (8 = 1, 2) (3.30) 

J U 

B=[I 2 -C 2 G,(4 1 ), & 2) )]- 1 C 2 [exp {a,?> {t- 1 )}^^), &«) + €,(#>, £>) exp {•^(r>)}*S?(4 1) , &»)] 

Ci(4 1 ),z)=ar 1 [exp{a 1 r^(t- 1 )}-I 1 ]a 2 ;C2 = a 4 - 1 [exp{a4? ! J(t- I )}-l2]. 

//ere Ii designates the kX k unit matrix, andlo designates the(n — /c)X(n— k) unit matrix, k being the 
number of elements in Ni (Section 3.2). All integrals and variations on the right of (3.30) are taken 
along SP. The bound on the right of (3.30) is valid provided that every eigenvalue of the matrix 
C-iCi^j, , £ k 2) ) is less than 1 in magnitude. Moreover, an exactly similar result holds for e^ k (z). 
Let us now obtain a norm bound. We define 

.//,>, v)=\\u^\\n, v (t-"')+ ti ^ > l^ nvy>(« 2 , u)p\At- m -') (3.3D 

where t\ and t 2 are any two points on $P in the order ($\ t\, tz, £ ( A 2) for 5= 1 and in the order @£\ 
t>i, t\, £ ( A 2) for s = 2. Similarly, with the notation of (3.27) we put 

a,= t S ^\\a k (t)\\ (k=l, 2, 3, 4). (3.32) 

Thus we obtain an inequality similar to (3.29); using Lemma 6 and the inequality (32) of the 
appendix we then obtain 

THEOREM 3.4: Let a k (k=l, 2, 3, 4) be defined by (3.32), let 0£>(k=l, 2) be defined by (3.31), 
and let 

f(z)=r c . 2 (t-') 

(3.33) 

F(z) = . a f : ' v> [e< a ' +a <> f < z > - 1 - (a, + a 4 )f(z)] 
\3.\ ~r 3.4) 

where the variations are taken along &. 7/F(£ ( k 2) ) < 1 then the norm of the vector €^ k (z) (equation 
(3.29)) is bounded by the quantity 

(1 -F(£< k »))-'{e a ' f < z >[F(z)+ 1 -F(a]W{L", *) 

+ a _f 3 [1 - e -< a ' +a ^ f < z >] [ e «" + 1 "« 3l »; - e< aj+a< > f < z )]^)(z, & 2 >) 

+ -^J- e ■*>+-4«'> [1 - e -< a ^ a <>™] ^>(£ k >>, £»). (3.34) 
ai ~r a4 

Moreover, an exactly similar result holds for e { ^ k (z). 

It is noteworthy that the conditions of Theorems 3.3 and 3.4 can always be satisfied by taking 
contours 8P sufficiently far away from the origin. 
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4. Appendix. Some Inequalities for Volterra and Fredholm Integral Equations 

4.1. Bellman's Lemma in Matrix Form 

We begin with an extension of Bellman's Lemma [5] to matrices. 

Lemma 1: Let <p, G and H be k X 1 vectors and E and F be k X k matrices of continuous non- 
negative functions of a variable cr in an interval a ^ cr ^ b such that F(o") E(<j) and I F(s)E(s)ds 
commute. If 10 ^ 

\<p(cr) ^ E(<r) J [F(s)p(s)+G(s)3ds + H(<r) (a^o-^b) (1) 

then the following inequality also holds in u^cr^b: 



tp(&) ^ H(o-) + E(o-) J exp |j F(T)E(T)dr| {G<s) + F(s)H(s)}ds. 



(2) 



The proof of Lemma 1 is a straightforward extension of the one-dimensional result in [5, 
page 135]. 

Lemma 2: With the conditions of Lemma 7, let E -1 (cr) and ^ a>b (E _, H) exist for a in (a, b). 
Then 



,<p(o-) ^ E(o-) exp [ J F(r)E(r)dr] |e '(a)H(a) + / a , , r (E ] H)+ J G{s)ds 



(3) 



To prove this result we integrate the right of (2) by parts. 

If E- ! (*)H(s) increases monotonically over (a,6), then / „, (7 (E- 1 H) = E- , ((r)H(o-)-E- , («)H(a); 
in this case the inequality (3) may be simplified. If, in addition, E(s) is a scalar quantity then (3) 
reduces to 



J> 



<p(&) ^ exp F(r)E(T)dr\ Hicr) + E(cr) G(s)ds\. (4) 





4.2. Bellman's Lemma for 


Fredholm Vector Integral Equations 






The following lemma e 


xtends the above 


results to Fredholm 


integral eq 


uations 






Lemma 3: When the variables cr, t, and 


t lie in the interval (a 


, b), let X(cr 


r)(kx 


k),rfr)(k 


XI), 


and 0{r) (k X 


1) be matrices 


of non-negative continuous functions 


such that 












X(o-, r) X(t 


, t) ^ X(o-, t) X(r, t) 








(5) 


and 




















<p(o) ^ J > X(o-, TV(T)dr+ 0((t). 








(6) 


Suppose that each eigenvalue of the matrix 
















F = 


f b X(r, r)dr 

Ja 








(7) 


is less than 1 


in magnitude. 


Then 














<p(o) 


^6(o)+\ X(cr 


THI-D-'flMdXa 


^ o- ^ b). 






(8) 


10 If (a, b) is not c 


)mpact we assume that t 


tie integral on the right of (1) 


exists and that each element of the 
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resulting vector is i 


iniformlj bo 


unded. 





PROOF: From the conditions of the Lemma there exists a unique kXl vector ^(cr) of non- 
negative continuous functions such that we can replace the inequality (6) by the equality 



<p(o) = I X(o-, T)(f(r)dT + 0(a) - V(a), 

J a 



(9) 



The equation (9) is a Fredholm integral equation in <p which we can solve by successive approxi- 
mations. Defining a sequence {<p v }(v = Q, 1, . . .) by <p = O, 



(10) 



<p v +i(<r)= f b X(o-, T)<pAT)dT+d(a)-V((r), v = 0, 1, . . . , 

Ja 

and using (5), we easily establish by induction that 

O ^ <p v+i (o)-<pJ(a) ^ [ b X(o% T)F"-H0(T)-V(T))dT 9 p= 1, 2, . . . . (11) 

If the eigenvalues of F are all less than 1 in magnitude, we have 

(I-F)-»=jF" (12) 

and clearly the sum of the power series on the right is a matrix of non-negative elements. On 
summing the inequalities (11), we obtain (8) with replaced by — ^, since it is well known that 
this sum bounds the true solution of (9). Decreasing the magnitude of the elements of W on the 
right of this sum increases the elements of (p. The inequality (8) now follows. 

The proof of the following Lemma is similar to that of the above Lemma, and is omitted. 
Lemma 4: With the conditions of Lemma 3 let (5) be replaced by 

X(o-,t)^G(ct)H(t) (5') 

where G and H are k X k matrices. Let each eigenvalue of 

F = f H(t)G(t)cLt (7') 

Ja 

be less than 1 in magnitude. Then 

<p(o) ^ 0(a) + I b G(o-)(I -F)- 1 H(r)^(r)dT (a ^ cr ^ b). (8') 

Ja 

4.3. Extensions to a Simultaneous Pair of Volterra Vector Integral Equations 

Here we consider the pair of inequalities 



<pi(o) ^ I {[oupi(s)^fi<pis)]w(s) + 1 (s)}ds 

J a 

<piar) £ f { [y<? 1(5) + 8<p-As) Ms) + 6. 2 (s)}ds 



(13) 



where </?i (kXI), <p 2 ((n — k) X 1), d\ (kXI), 6 2 ((ft -r #c) X 1), co (1X1) are matrices of nonnegative 
continuous functions on (a, b), and n ^ 2, < k < n. Also, a, j3 9 y, and 8 are constant matrices 
of orders kXk, kX(ti — k), (n — k)Xk and (n — k) X (n — k) respectively, whose elements are 
nonnegative. 
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On applying Lemma 1 to each of the equations (13), we obtain 

Ho) SJ ™ s) -* A<r) {y<pi(s)a)(s) + Ods)}ds 
where 

f(s)=[\{r)dT. 
Substituting the second of (14) for <£> 2 in the first of (14), we get 

^ l(0_) ~ P I" eaf(a) ~ af(s) \fi e df{T) - 8f{s) (Wi(t)co(t) + 2 (r) }drco(s) + 0i(s)]ds. 



(14) 



(15) 



(16) 



The repeated integral (16) may be interpreted as a double integral over the shaded region indicated 
in figure 3. Splitting this region into a triangle and a rectangle and interchanging the order of 
integration, we obtain 

tpi{&) ^ 3 (tr) + 1 J e*>W-*W pe W-%w ly^^r) + 0>(t)}a)(s)dsdr 

A h \\ af ^- af(s) p e * fiT) -* fis) {y<pi(T)€o(T) + d.2(T)}a>(s)dsdT, (17) 



where 



+ 



fe(cr)= I e ««»>-*« ^( s )d5. 



(18) 



Figure 3. 
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We evaluate the inner integrals in (17) to find that 

Wo-) g 6» 3 (o-)+ I X(cr, T){y<p,(TMT) + 0-2(T)}dT (19) 

where X(cr, r) satisfies 

f e a/(CT) [pe 8 ^ (T) - e " <^ (t) 0] , r ^ cr 
aX(cr,T) + X(cr,T)8=| (20) 

The equations (20) have unique solutions for X(cr, r) provided that /Xj-h Vj^ (£=?= 1, 2, . . . , 
k; ,/== 1,2, . . . , n — k), where /x* and i>j are the eigenvalues of a and 5, respectively. Moreover, 
this condition can always be satisfied by an arbitrarily small change in the elements of a and 5, if 
necessary. Nevertheless an exact explicit expression for the right hand side of (20) is generally 
not feasible, and in practice we make a slight sacrifice of sharpness of error bound in order to 
achieve a matrix bound which is simpler to evaluate. Several simplifications are possible; one 
of these is as follows: n 

e a/(cr)- a f(s) fie^M-m*) (jj(s)ds^ J e«/<°-)-«/< 5 ) (5e 8 f^ ais)ds 

J a J a 

= a \e afM - e **•»-««•> ]p e 8/w Si a~\e ^ - 1 K ]j3e S/,T) . (21) 



(22) 



Similarly 

f <r e «w-4»0e«M-«u M{s )ds g or* [ e " M) - I K ]/8e m - G(o-)H(t), 

where 

G(a-) = «-'[e^-I K ], H(r) = /3e s ^. (23) 

Combining (22) and (23), we obtain 

W<t) 5= 03(0-) + f " G(o-)H(t){ W .(t)w(t) + Wt)} *. (24) 

J a 

We apply Lemma 4 to this inequality to obtain: 

Lemma 5: Let <p u <p 2 , «, j3, y, 8, co, 0, an^ 2 6e defined as in (13) ana 7 the functions G and H 
as in (23). //Jne inequalities (13) Zio/a 7 ana 7 i/eacn eigenvalue of the matrix 



'=[ b H(r) y G( 

J a 



r)w(r)dr (25) 



is less than 1 m magnitude, then 

<pi(o) ^ 6(<t)+ f G(o")(I-F)- 1 H(s) 7 ^(s)w(s)ds, 

Ja 



(26) 



"Since a _1 [exp(a) — I] = V ,,,,,, the possibility of a being singular is not excluded. 
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where 



0((t) = OA(t) + r G((T)H(r)0-AT)dr (27) 



and Ox is defined by (18). Moreover an exactly similar result holds for (p^cr). 

In the case when n = 2 and k=\ the equations (20) can be solved explicitly. We then have 
(19) with 



r e 8/(r) _ e -a/(r)"| T = O" 



X(a, r)= (28) 

It is easily verified that for arbitrary points cr, r and £ on («,&), X(o\ t)X(t, /) = X(cr, t)X(r, r); i.e., 

the conditions of Lemma 3 are satisfied. On applying Lemma 3 to the inequality (19) we obtain 

Lemma 6: Let ip u <^ 2 , «, /3, y, 5, co, 0i and 2 be defined as in (13), wi7/i n = 2, k= 1, arcd 

F(cr)= ^ y [e ( " +fi)f(<r) -l-(a + 8)f(cr)], (29) 

where f(s) is defined by (15). If the inequalities (13) Ao/rf a/id (f F(b) < 1, j/iera 

?i(o) ^ ^°") + 1 _F (b) [" X(cr ^ T)[ya)(T)^(r) + ^ 2 (t) ]dr (a ^ <r ^ b) (30) 



wAere X(cr, r) is defined by (28) and 6 :i is defined by (18). 

On expanding the right hand side of (30) and using the inequality X(cr„ t)X{t, t) ^ X(cr, f)X(r, t) 
once more, we get 



R(> CtfitT) 



l-F(/;)(a + 8) 



| fr j^^ [F(<r) - F(5) + 1 - F(6)]e- a/w » 1 ( 5 )+ [e^ s) - e"^]^)} efc 



(31) 



For ease of evaluation we make some over-estimates in (31) to obtain 



!«*» [F(o-) 4- 1 - F(6)] J 9i(t)* 



2 

L-F(A) 



+ ^_ C *to+W) n_ e -<«+«Wo)] P ^(rV/rl- (32) 

a + o J (l J 

In conclusion, we observe that if we replace all inequalities in Lemmas 1,2,3, and 4 by equali- 
ties, then the results of these lemmas are exact; hence the inequalities are sharp. However, since 
X(<r, r) < G(cr)H(r) almost everywhere in a =i cr, r = b the result of Lemma 5 is not sharp. The 
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result of Lemma 6 also is not sharp because, although X(cr, t)X(t, f) = X(cr, J)X(t, t) for <x, r, £ in 
the order o% r, £ or t, r, cr on (a,6), X(o% t)X(t, < X(cr, t)X{r, r) when this order of o% t and * is 
violated. 
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